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ABSTRACT 

The snag of common linear filtered-G Adaptive 
Inverse Control (AIC) method is modify to hold with 
the characteristic of non-linear object with time delay 
and corresponding filtered G-adaptive Algorithm 
based on Real-Time Recurrent Learning (RTRL) is 
presented to identify the parameters and design the 
controller. The simulation result on a non-linear ship 
model of “The R.O.V. Zeefakkel” and Adaptive PID 
control keeps the same dynamic response 
performance, and also mathematical model for ship 
Course keeping be discussed. A neural network 
adaptive filter is introduced for the removal of 
impulse noise in digital images. 

Keyword: Filtered- E Adaptive Inverse Control, Ship 
course keeping, impulse noise detection and removal. 

INTRODUCTION 

Adaptive Inverse Control (AIC) is novel approaches 
which can make a plant track the input command 
signal with a controller whose transfer function 
approximates the inverse of plant transfer function. 
Compared with traditional methods, AIC can achieve 
specified dynamic responses more easily and has 
better ability of disturbance rejection. Simulation 
comparison with previous scheme and adaptive PID 
control are performed on the non-linear ship 
maneuvering to test the effect of improved algorithm. 
The impulse noise detection and removal of impulse 
noise while preserving the integrity of an image 
expressed mathematically be discussed. 

The Improved Filtered-G Adaptive Inverse 
Control 

The Single Input Single Output (SISO) discrete-time 
nonlinear system to be considered is described by n 
dimension state equations as below 


x(k + 1) = f[x{k),u(Ji)], f (0,0) = 0 

y(k) = h[x(k)] 1 } 

Where /(.) is a non-linear function, u(k) and y(/t) 
respect output and input variables, respectively. The 
principles of inverse control can be extended to deal 
with non-linear system, through non-linear system has 
no strict inverse model. Whereas a linear system 
possesses a unique inverse, non-linear systems have 
only local inverses, valid in a bounded region of the 
signal space. The non-linear system is supposed to has 
Bounded-Input Bounded Output(BIBO) stability, i.e., 
there are existing constants C ± and C 2 , and for the 
input u(k) G s u ,s u = {u\\u(k) |< C 1; V/r > k 0 }, the 
system output y(/c) satisfies: |y(k)|< C 2 . If its 
linearized system is observable at balance point 
x = 0, the non-linear system described by (1) has 
local NARMA (Non-linear Auto-Regressive Average) 
model as follows 

y(k + d) = F[y(/t ),... 
y{k — Uy), u(Ji),..., u{k — n u )] (2) 

Where n u and n y are the orders of input time series 
and output time series respectively: n u < n y , d is 
system delay which can ensure the existence of 
system inverse model, and F(.) is certain non-linear 
function. If F(.) is mono tonic change with (k) , 
system (l)is invertible in the input spectrum S u . If the 
input signal is y m (k + 1) at time(/c + 1), the 
controller can be expressed as 

u(k) = g[y m (k + 1 ),y(k), . ,y{k - n y + 

1, uk~ 1, . uk~ nu+1 (3) 

Function F( .) and g( . ) are both unique. Under the 
dreive of u(k), the plant output will track the input 
signal of the inverse controller. 
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Identification of Delay Time Complexity 

The output of a non-linear system with time-lag 
characteristic can be described as below form 
y(fc) = f(x(k — d ),..., x(k — m — d ), 
y(k- 1),.y(/c — n)) (4) 

Where x(k) and y(/c) represent the input and output 
of system at time k, m and n are the input and output 
order, d is system delay time, and /( .) is the transfer 
function of this non-linear system. 

If we use neural network to model such a non-linear 
plant, the training effectiveness will be greatly 
influenced by the samples of input vector collected in 
different range. Therefore, such network can be used 
as delay-time identifier. Here a non-linear auto¬ 
regressive with exogenous inputs ( NARX) network 
is used, which can be described as 

x(k) = [x(k),... .,x{k — a), 
y(k-l),....y(k-b)] (5) 

Where a and b denotes the length of tapped delay 
lines of input and output feedback, respectively. The 
plant input/output data pairs collected in certain time 
period are used to train the NARX network, and many 
research result have demonstrated that the sample 
range of u(k) changing from not including the first 
delayed input signal a = d — 1 to including it, 
i.e.,a = d — 1 -» a = d, will cause the training 
effectiveness having a great improvement, and hence 
the delay time d of controlled plant can be confirmed. 

Improved Filtered-£ Algorithm 

The NARX filter is also used to the controller 
modeling. For including self-feedback, it must be 
adapted using a method such as Real-Time Recurrent 
Learning (RTRL) or Back-propogation Through Time 
(BPTT). Considering the real-time requirement of 
practical application, we propose the filtered-^ 
algorithm based on RTRL. To model C, RTRL usually 
needs to obtain the error £* between C and the ideal 
controller outputs, and then C parameters can be 
adjusted via the principle of gradient descent to 
minimize the objective function ) c = (£*) 2 /2 , but £* 
is unknown in practical AIC implementation. In 
improved filtered-£ algorithm, £* is replaced by 
system error £ which is filtered by 1 to adapt C 
parameters. The update algorithm can be expressed as 
follows 

W c (k + 1) = W c (k ) 
+H(l-r 1 )(-y c (k)) + r 1 W c (k) (6) 


Where /i and rj denote the step size and momentum 
co-efficient respectively, and V c (k) represents the 
gradient of J c with respect to controller weights W c : 

V(k ) = V£ = 10^31 

c v y wt/17 n 


= £(k) 


dW c 2 dW c 
d(u(k)-u(k- A) 


dW r 


= -s'(k) 


du(k- A) 
dW r 


( 7 ) 


Where ~ is the A —delayed partial derivative of 

C output u(k) with respect to its weights w c at time k. 
It is updated by the chain rule: 
du(k ) 


dW f 


C 

du{k ) 
dW r 


u 

I 

n=0 


du(k ) dr(k — ri) 
dr(k — n) dW c 


+ 


I 

771 — 1 


du(k) du(k — m) 


du(k. — m) dW c 


_ du(k) , v->i) du(k ) du(k-m ) 

- V hm=\ —— (8) 


dWr 


- du(k-m ) dW c 


The terms d ^ k \ and are the direct 

dW c ’ dr(k-n) dQ(k-m) 

effect of a change in the weights, previous inputs and 
outputs on the current output respectively, and these 
partial derivatives may be computed directly using the 

back-propagation algorithm. dw — is simply a 


du(k) 


previously-calculated and stored value of —— . For i 


= 0 , - 1,-2 . 


du(i ) 

dWr 


dW c 

are set to zero. We also note 


dv(jc 71 ) a 

that dw —, the partial derivative of the input of C 
with respect to its weights, is zero for allk. 

SIMULATION RESEARCH 

We apply the improved Filtered-£ AIC method 
proposed in this paper to the control of non-linear ship 
maneuvering. As a typical non-linear system with 
time-lag characteristic, ship can be described by the 
non-linear model. 

Tip + ip + aip 3 = KS + a) (9) 

Where xp is the yaw angle, 6 is the control rudder 
angle, T,K and a are the model parameter, and co is 
the external disturbance. 


The ship model used in the simulation study is 
adapted from Van Amerongen (1982). In this 
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reference , “The R.O.V. Zeefakkel” is described by 
the following set of parameters and the ship ‘R.O.V. 
ZeefakkeT of 45 m length can be represented by the 
non-linear mathematical model and the model 
parameters for a ship speed of 10 knots be Ship length 
L = 45m, service speed v = 5.0m/s, gain constant 
k = 0.5s -1 , time constant T = 31s, and Norrbin co¬ 
efficients n x = 1.0, n 3 = 0.4(s) 2 . In order to evaluate 
the performance of the improved filtered-£-AIC 
scheme, conventional filtered-^ AIC and adaptive PID 
control are also performed. In the simulation, the 
sampling period 

T = 0.1s. 

(a)Firstly the ship +10° course-changing test is 
performed without external disturbance, namely 
a) = 0. At k = 300s, the ship time constant constant 
T changes from 31 s to 45 s. Clearly the conditional 
filtered-^ AIC has much better dynamic performances 
than that of the adaptive PID control, but there 
existing static error in it for the influence of the 
change of plant delay time. In contrast, the improved 
scheme avoids static error and at the mean time, it has 
kept the same dynamic performances as conventional 
filtered-^ AIC method. 

(a) The second simulation to be examined is the 
disturbance rejection ability. The external disturbance 
is chosen as 

a) = 4.5X + 3.44 °H 2 (10) 

Where H 1 and H 2 are normally distributed random 
numbers with mean 0, varience a 2 = 1. The 
comparison of the response curves between improved 
method and conventional ones. It can be clearly seen 
that the improved scheme has better disturbance 
rejection ability. 

Ship Course-Keeping Mathematical Model 

Equation (1) gives linear first-Order Nomoto model 
which is used for ship course-keeping 

yP+±*P = ^s, (1) 

! 0 ! 0 

Among these, ip is the ship heading angle, S is the 
rudder angle, K 0 is the steering gain index, T 0 is the 
ship following index. It presents a Nomoto model 
based responding type of non-linear mathematical 
model, that is using a non-linear item (K 0 /T 0 )H(ip) to 
replace ip/T 0 , and 

H(xP)=xIj + Pi, 3 , ( 2 ) 


Where a, /? represent the scale co-efficient of the first 
and the third power of the rate of turning i p . Specific 
parameter value depends on the factors like ship type, 
stowage, ship speed, etc. Taking xp as ship heading 
angle, then ip = r; meanwhile, taking the frequently 
influence from uncertain interference such as wind, 
wave and current into consideration, the non-linear 
dynamic equation of ship course-keeping system 
could be written as: 

I i p = r 

f = -j-(ar + pr 3 )+ ^-8 + A, ^ 

! 0 1 0 

Where A is the uncertain interference. According to 
engineering practice, A is usually bounded 
disturbance, assuming there is an unknown constant 
could meet the equation (1) on the basis of hypothesis. 

||A||oo< P , p is a unknown positive constant. 

( 4 ) 

Impulse Noise Detection and Removal 

When an image is coded and transmitted over a noisy 
channel or degraded by electrical sensor noise, as in a 
vidiction TV camera, degradation appear as “salt 
and-paper” noise (i.e., positive and negative 
impulses). Two models have been proposed to 
describe impulsive noise (Justusson, 1981). The first 
model assumes fixed values for all the impulses, and 
the second model allows the impulses to follow a 
random distribution. In this work, only the first model 
is dealt with. The model is mathematically expressed 
as 

{ d p with probability P p 
d n with probability P n 
fi with probability 1 — (P p + P n ) 

Where f t denotes the uncorruped pixel values, d p 
denotes the fixed value of the positive impulses, d n 
denotes the fixed value of the negative impulses, and 
x t denotes the pixel values of the degraded images. 

Noisy pixels can be characterized by their local 
statistical properties. To extract features from local 
statistics, a window is used to pass through the entire 
noise degraded image. Theoretically, the size of the 
window is arbitrary. However, a 3x3 window is good 
enough in most applications. From the small window, 
useful local features are obtained, such as mean, 
median, variance, range, and extreme range. Noise 
detection is then treated as a pattern classification 
problem. Assume that M features are measured from 
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each input pattern. Each set of M features is 
considered as a vector in the M-dimensional feature 
space. Pattern classification is accomplished by a self- 
organising neural network which is capable of cluster 
analysis for image data. 

The neural network proposed in this paper is a 
simplified version of ART I (Carpenter & Gross berg, 
1987). However, it is capable of grey - level 
classification. It has two layers: the input layer and the 
output layer. These layers are connected by feed 
forward paths. The nodes of output layer are lateral 
inhibitive to one another. They are designed as a 
competitive network capable of choosing the winning 
node. The set of weights u £; - for any specific node 
i constitute the prototype pattern for cluster!. The 
clustering algorithm starts from isolated patterns and 
coalesces the nearest patterns or groups according to a 
threshold from the bottom up to form hierarchies. The 
essential point of this algorithm is to build up the 
clusters using the Euclidean distance measure 
between the input s and the weightsiq. Assume a 
sequence of pattern samples 

s = (s 1 ,s 2 ,....s M ) e R m 
and a set of weight vectors 

Ui = (u il ,u i2l ...,u iM ) E R m 
It is well known that similarities exist between ART 
and Kohonen’s self-organising map. 


CONCLUSION 

In this paper, an improved nonlinear filtered-£ method 
is introduced and adaptive PID control are performed 
on the control of non-linear ship maneuvering and 
also mathematical model for ship Course-Keeping and 
impulsive noise be expressed. 

REFERENCE 

1. Widrow B, Walach E. Adaptive inverse control. 
Englewood cliffs, N J: Prentice-Hall press, 1996. 

2. Zhang Y H, Liu S X. An adaptive automatic ship 
steering system. Shipbuilding of china, 1991. 

3. Nieminen, A., Heinonen, P., & Neuvo, Y. (1987). 
A new class of detail-preserving fdters for image 
processing. IEEE Transactions pattern Analysis 
Machine Intelligence. 

4. Zhang Xianku, Lu Xiaofei, Guo Chen, Yang 
Yansheng, Robust Neural Network Control for 
Ship’s Course-Keeping, Journal of Ship 
Mechanics, Vol. 10, Issue 5, 2006. 

5. Du Gaang, Zhan Xingqun, Zhong Shan, Zhang 
Weiming, Adaptive inverse control of non-linear 
ship manueuvering based on Neural Networks, 
Measuring and Testing Technique, Vol. 24 Issue 
4, 2005. 


@ IJTSRD | Available Online @ www.ijtsrd.com | Volume-2 | Issue-5 | Jul-Aug 2018 


Page: 1994 



